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Free standing InP quantum dots have previously been theoretically and experimentally shown to have a direct 
band gap across a large range of experimentally accessible sizes. We demonstrate that when these dots are 
embedded coherently within a GaP barrier material, the effects of quantum confinement in conjunction with 
coherent strain suggest there will be a critical diameter of dot (~60A ), above which the dot is direct, type I, and 
below which it is indirect, type II. However, the strain in the system acts to produce another conduction state 
with an even lower energy, in which electrons are localized in small pockets at the interface between the InP dot 
and the GaP barrier. Since this conduction state is GaP Xi c -derived and the highest occupied valence state is 
InP, F-derived, the fundamental transition is predicted to be indirect in both real and reciprocal space ("type II") 
for all dot sizes. This effect is peculiar to the strained dot, and is absent in the free-standing dot. 
PACS:73.20.-r, 73.20.Dx, 85.30.Vw 



There are two leading techniques for fabricating InP 
quantum dots; (i) colloidal growth, producing un- 
strained, chemically passivated dots 0] and (ii) dots 
grown by the Stranski-Krastanov mechanism that are 
embedded within a semiconductor barrier such as GalnP 
. In a recent study of the electronic structure of col- 
loidally grown InP dots |Q it was found that these dots 
have a direct band gap at the T-point of the Brillouin 
zone for all experimentally accessible sizes. In contrast, 
it has previously been shown that for free-standing GaAs 
dots, the band gaps can undergo a transition from direct 
to indirect as a result of quantum size effects §. The 
effects of quantum confinement are that as one decreases 
the size of the dot, all of the conduction levels are pushed 
up in energy at a rate reflecting approximately the in- 
verse of the electron effective mass. Since the Ti c masses 
are generally lighter than X\ c masses ||, reduced sizes 
can transform a direct dot into an indirect dot if the ini- 
tial Yic-Xic separation in the bulk is not too large. It 
has been predicted that in GaAs, where the bulk Ti c - 
X\ c separation is only 0.55 eV §, a free-standing, zero- 
pressure dot will become indirect at a dot size of 40A , 
whereas an AlAs-embedded GaAs dot will become indi- 
rect at 80A 0. Since, however, in InP the Ti c -Xi c bulk 
separation is large (0.94 eV ||), calculations ||] have pre- 
dicted that this separation is not overcome by quantum- 
size effects, and free-standing dots will remain direct at 
all sizes at zero pressure. 

In this paper we are interested in investigating whether 
InP quantum dots embedded within a GaP matrix ex- 
hibit a direct-to-indirect transition as a function of 
size and thereby exhibit significantly different electronic 
structure to free standing InP quantum dots that are 
direct for all sizes. We find indeed that for spherical 



InP dots smaller than 60A in diameter, the X± c conduc- 
tion band of the unstrained GaP barrier is lower than 
the Tic-state of the InP. This is analogous to a "type 
II in real-space and in reciprocal space" state familiar 
II in Al As/ GaAs nanostructures. Surprisingly, however, 
we find that under the influence of the hydrostatic and 
biaxial strain present at the GaP/InP interface, a new 
conduction state emerges that is lower in energy than 
both the unstrained bulk GaP Ai c -state and the InP Tic- 
state. This qualitatively new type of state is localized at 
the interface of the dot and its barrier, and is indirect 
(A-like). Hence, when coherency exists between InP and 
GaP, we predict that photo-excited electrons will be lo- 
calized in this state, giving rise to an unusual dependence 
of the band gap on size. This effect is peculiar to co- 
herently strained systems, and is absent in free-standing 
(colloidal) dots. 

Expected trends based on band offsets: Before presenting 
our calculated results for GaP-embedded InP dots, we 
discuss the basic expectations regarding the nature of the 
confined states. Figure |l|a shows our fitted |)| unstrained 
( "natural" ) valence and conduction band offsets between 
bulk InP and bulk GaP. The T states are shown as heavy 
solid lines and the X states as thin solid lines. We see 
that InP can act as a "well" both for the conduction band 
Tic electrons and the valence band Ti 5 „ holes (a "type I" 
offset). The confined levels are denoted schematically by 
dashed lines. As the InP dot becomes smaller, quantum 
confinement causes the confined InP valence levels to be 
pushed down in energy and the confined conduction levels 
to be pushed up (see arrows). This causes the Tic-level 
to be pushed up in energy with respect to the Xi c -level 
in the GaP barrier causing a direct-to-indirect transition 
(see below). 
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FIG. 1. (a) Unstrained ("natural") band offsets (in eV) 
between bulk GaP and InP. Solid lines indicate bulk band 
edges and dashed lines indicate quantum confined levels. Ar- 
rows show the energy change due to confinement. P-derived 
states are shown with thick lines and X-derived states with 
thin lines, (b) Strain modified band edges, _E n k[e(R)], plotted 
along [100] through the center of the InP dot with diameter 
13lA and dot-dot separation of 109A. The lowest(highest) 
conduction(valence) band is shown at each position R. The 
* denotes the position at which the lowest conduction state 
is localized. 

In reality, the large (7%) atomic size-mismatch be- 
tween GaP and InP will cause atoms to be displaced 
off their ideal zincblende positions. This will alter the 
effective band offsets of Fig Qa and thus the anticipated 
confinements. To calculate the ensuing strain we place 
an InP sphere of radius R inside a large GaP cube, and 
relax all the atomic positions to their minimum strain 
energy values, using the Valence Force Field (VFF) elas- 
tic energy functional |f[o|| . We chose to fix the external 
dimension of the GaP cube during the relaxation as this 
most closely resembles the experimental situation where 
InP dots are grown (|] on a fixed GaP substrate and the 
dot-dot separation is large enough to remove any dot-dot 
interactions. However, in the systems studied here, the 
large barrier sizes create such a high GaPTnP ratio (~ 
30:1) that any external relaxation would be minimal in 
any case. The resulting strain exhibits nontrivial hydro- 
static and biaxial components. Our quantum mechani- 
cal calculation of the energy levels of the dot (see below) 
will include the effect of such a strain profile. However, 
in order to understand these results, we hrst consider a 



simpler case, namely we calculate the band edge states of 
bulk InP and bulk GaP subject to the local strain, e(R), 
experienced by the GaP embedded InP dot at position 
R. To do this we discretize the GaP/InP nanostruc- 
ture into "cells" with position vector R and then per- 
form ^40 bulk band structure calculations of InP and 
GaP, using the empirical pseudopotential method | pl| , 
thus obtaining the bulk eigenvalues i?„fe[e(R)] for band 
n at wavevector k within each cell. Each bulk calcu- 
lation, i?„fc[e(R)], uses the In-P or Ga-P bond geome- 
try within that cell. The resulting strain-modified band 
edge states are shown in Fig. |j]b. Compared with the 
unstrained offsets (Fig. |l|a), we see that the GaP X\ c - 
band edge that is flat in the absence of strain (Fig. |l]a) 
is now transformed into an attractive trough (indicated 
by * in Fig. |l|b), capable of localizing electrons. The 
formation of this trough is initially surprising as the de- 
formation potential at the X lc -point is positive and one 
might therefore expect the hydrostatic expansion of the 
GaP at the interface with the InP dot to drive the X\ c - 
state up in energy. However, the above bulk calculations 
show that it is the biaxial strain present at this inter- 
face which is the dominant term, and this is capable of 
forming the electron troughs. The atomistic strain has 
therefore profoundly modified the nature of the confined 
electron states from delocalized to localized. It is im- 
portant to emphasize that conventional |l2| ] calculations 
of strain modified conduction band offsets include only 
the hydrostatic (no biaxial) term and only the Tic (no 
X\ c ) conduction band, and would therefore miss the im- 
portant changes in the conduction band edges between 
Figs. [j]a and 0b which our calculations show are due to 
the effect of biaxial deformation on the X\ c state. 
Results of calculations on dots: To calculate the energy 
levels of GaP embedded InP dots, we again place an InP 
dot of radius R, surrounded by sufficiently thick GaP 
barrier in a "supercell", repeated periodically to create 
a lattice of dots. Having created (artificial) translational 
periodicity, band theoretical methods can then be ap- 
plied to study the electronic properties. The limit of an 
isolated dot is achieved by increasing the thickness of the 
GaP barrier. The calculations for both the bulk bands 
and the quantum dot levels are based on the atomistic 
Hamiltonian 



H 



Rn 



The total potential is constructed from screened atomic 
pseudopotentials, v a , where a represents Ga, In and P, 
and R Qn are the relaxed atomic positions. The pseu- 
dopotentials, v a , have been fitted Q to the experimental 
band gaps, deformation potentials and effective masses. 
We use the analytic form of the pseudopotential de- 
scribed in Ref. [[J, which was designed to build in the 
effects of strain experienced by each atom in lattice mis- 
matched systems. 
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The supercells studied in this paper contain up to one 
million atoms, which is too large for the Hamiltonian in 
Eq.(|l]) to be solved by direct diagonalization. We thus 
use the "Folded Spectrum Method" (FSM) in 
which one solves for the eigenstates of the equation 



H - e 



ref 



ij>i = (e- e ref ) 2 ^i 



(2) 



where e re r is a reference energy, and the wavefunctions, 
ipi, are expanded in a plane wave basis. By placing e re f 
within the gap, and close to the valence band maximum 
or conduction band minimum, one is then able to obtain 
the top few valence states or the bottom few conduction 
states respectively. Using this approach the computa- 
tional cost scales as MNlog(N), where N is the number 
of desired electronic states and M is the number of plane 
wave basis functions (M « 20 million in the largest sys- 
tem studied here). The simulations in this paper were 
performed using a parallel code on the Cray T3E900 on 
up to 256 processors where the ellipsoid of g- vectors is di- 
vided over the processors in a similar way to the method 
used by Clarke et al Using this data distribution 

and fast parallel FFT's, almost linear speed up with the 
number of processors can be obtained for the large sys- 
tems studied. 

To facilitate comparison with the spherical free stand- 
ing InP quantum dots studied in Refs. [Q, we constructed 
a series of supercells containing spherical InP quantum 
dots with diameters of 44, 87, 131 and 174 A. Each 
dot was surrounded with sufficient GaP barrier mate- 
rial to produce dot-dot separations of 109, 152, 196 and 
239 A respectively. The calculated energies of the high- 
est occupied valence states and lowest empty conduction 
states are shown in Fig. ||. The left hand side of Fig. || il- 
lustrates the corresponding wavefunctions squared of the 
131 A dot. We see that the highest energy valence 
wavefunction is localized within the InP dot (Fig. ||c), 
whereas the lowest conduction wavefunction is localized 
in pockets at the {001} facets of the interface between 
the InP dot and the GaP barrier (Fig. ||b). The energy 
of this interfacial state is considerably lower than the un- 
strained bulk GaP Xi c -state (solid line in Fig.^) for all 
the dots studied. To establish the identity of these wave- 
functions in terms of the parent GaP and InP bulk states, 
we project the dot wavefunctions, ipi, into the zincblende 
Brillouin zone using the method described in Ref. [115 
This mapping is shown on the right hand side of Fig 
for the k z = plane through the Brillouin zone. We 
see that the highest energy valence state is a T-derived 
state (Fig. ||c), while the lowest conduction state is X- 
derived (Fig. ||b). The calculated dipole transition ma- 
trix clement between these states is 5 orders of magnitude 
smaller than one would expect between a more typical 
pair of T-derived conduction and valence states, render- 
ing the transition forbidden. 
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FIG. 2. Energies of the near edge states of GaP-embedded 
InP dots with diameters 44, 87, 131 and 174 A and a dot-dot 
separation of 65A 

To investigate the effect of dot-dot separation, a series 
of similar calculations were performed on the quantum 
dot with a diameter of 87 A, where the dot-dot sepa- 
ration was steadily increased to values of 152, 196 and 
240 A. The energies of the A-derived, interface localized 
states were -4.15, -4.22 and -4.22 eV respectively. The 
invariance of this state with dot-dot separation strongly 
suggests that the strain induced interface localization ev- 
ident in Fig ||b is not due to too small a choice of super- 
cell, but an intrinsic effect present at such a interface. 
Given that the lowest-energy conduction state, as illus- 
trated in Fig. ||b is an Ai c -derived, interfacial state, one 
wonders where is the Tic-derived dot conduction state. 
It would not be practical to use the FSM to search for 
this state as the large size of the supercell folds many 
eigenstates to the T-point (at which the FSM calcula- 
tion is performed) between the lowest conduction state 
and the Tic-derived state. We instead use the Linear 
Combination of Bulk Bands (LCBB) Jl7|] method which 
specifically searches for states of a given symmetry (e.g. 
Tic-derived), whether or not they are the lowest energy. 
Here one first solves for a set of bulk Bloch wavefunc- 
tions, ipnki of the two materials, InP and GaP. Then the 
Hamiltonian from Eq.(|l]) is diagonalized within the basis 
of these wavefunctions. The LCBB method allows one 
to choose which bulk Bloch wavefunctions to include in 
the basis set. As we are only interested here in the Tic- 
derived energy levels, we include bulk Bloch wavefunc- 
tions from a radius of (in reciprocal space) around 
the T-point, where I is the supercell length. The resulting 
eigenstates were found to be converged with respect to 
the basis size at this radius. The wavefunction squared 
of the Tic-derived state for a system containing a dot 
with a diameter of 131 A is plotted in Fig. ||a. The en- 
ergies of the Tic-derived states in each of the four InP 
quantum dots are shown in Fig. ^. Figure. || shows that 
there is a critical dot diameter around 60 A below which 
the Tic-derived conduction state in the InP quantum dot 
is higher in energy than the bulk Ai c -state of the GaP 
barrier. This is a type I to type II transition. However, 
for all sizes of InP dot, the Ti c -derived state is higher in 
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energy than the Xi c -like interfacial state. 
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FIG. 3. Wavefunction squared (left) and momentum-space 
analysis (right) for the near edge states (see Fig. ^) of an InP 
dot with a diameter of 13lA and dot-dot separation of 109A. 
The left hand side shows each wavefunction squared in the 
(001) plane through the center of the InP dot. The right hand 
side shows the momentum-space projection of each wavefunc- 
tion in the kz = plane of the Brillouin zone. Wavefunctions 
(a) and (c) are T-derived, and (b) is X-derived. 

In conclusion, we have shown that: (i) the effects of 
quantum confinement and pressure raise the energy of 
the Tic-derived state in an InP quantum dot so that for 
dots smaller than 50A, this state is higher than the X\ c 
state of the unstrained GaP barrier. This transition is 
analogous to the AlAs-embedded GaAs dot, where the 
CBM moves from GaAs-Ti c to AlAs-Xi c as the GaAs 
size decreases. However, (ii) strain induces an even lower 
energy state, indirect in reciprocal space and localized in 
real space at the interface between the InP dot and the 
GaP barrier. Therefore, we predict that, even for large, 
spherical InP dots, as long as coherency is maintained, 
the effects of strain create a system with an indirect band 
gap, that is considerably reduced due to the low lying 
interfacial state. This is in direct contrast to the the 
behavior of free standing InP dots, which are direct over 
the large range of experimentally accessible sizes. 
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